Size and mass of Cooper pairs determined by low-energy /j,SH and PNR 
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The Pippard coherence length £o (the size of a Cooper pair) in an extreme type-I superconductor 
was determined directly through high-resolution measurement of the nonlocal electrodynamic effect 
combining low-energy muon spin rotation spectroscopy and polarized neutron reflectometry. The 
renormalization factor Z = m* cp /2m (m* p and m are the mass of the Cooper pair and the electron, 
respectively) resulting from the electron-phonon interaction, and the temperature dependent London 
penetration depth Al(T) were determined as well. A general expression linking £o, Z and Al(0) 
is introduced and experimentally verified. This expression allows one to determine experimentally 
the Pippard coherence length in any superconductor, independent of whether the electrodynamics 
is local or nonlocal, conventional or unconventional. 

PACS numbers: 



A core concept of superconductivity is Cooper pair- 
ing of electrons. Cooper pairs have a characteristic size 
£o, the Pippard coherence length, and an effective mass 
m* p = 2Zm, m being the Coulomb and band-structure 
effective mass of electron and Z the renormalization fac- 
tor, a measure of the electron-boson coupling strength [l[ . 
Only weakly dependent on temperature T, £o provides 
a reference length for the fundamental length scales, 
including in particular the London penetration depth 
Al(T), characterizing the decay of a penetrating mag- 
netic field. In units of £ , -V(T) determines whether 
the superconductor is described by local or nonlocal elec- 
trodynamics, and Al(0) determines whether it is type-I 
or type-II Q. To ensure consistency £o and Z should 
be measured simultaneously. However, in spite of their 
importance, £o an d Z have not yet been directly and si- 
multaneously measured in any superconductor. 

As shown by Pippard [|| , £o can be determined in non- 
local superconductors through measurement of the non- 
local electrodynamic effect [J] . Most superconductors are 
local and for those it was stated that £o cannot be mea- 
sured directly @. Here we report on a direct and si- 
multaneous determination of £oj Z and Al in a nonlocal 
superconductor (In) and answer the question how £o can 
be determined experimentally in any superconductor. 

When a superconductor is in the Meissner state, an 
external magnetic field Bq is completely screened due 
to a persistent current running in a thin surface layer 
over which the field decays to zero. The layer thick- 
ness is of the order of the "magnetic penetration depth" 
A = Bq 1 Jp 00 B(z)dz, z being the distance from the sur- 
face. If the size of the Cooper pairs is small (£o <C A), the 
relationship between the current density and the vector 



potential can be treated as local. Then the field decays as 
cxp(— z/Al), where Al = yj 'Ac 2 /Att is the London pen- 
etration depth (c is the speed of light). In the London 
theory the phcnomcnological parameter A is a function 
of the mass and the density of superconducting electrons, 
none of the two being well defined Q. 

The local approximation is applicable to superconduc- 
tors with Ginzburg-Landau parameter k > 1.6 @. For 
superconductors with smaller k the Cooper pair size is no 
longer negligible and the current density is determined by 
the vector potential averaged over a region of dimension 
£o- The nonlocality leads to deeper field penetration and 
to distortion of the B(z) shape. In the pure limit (elastic 
mean free path I 3> £o) B(z) is a function of both \l(T) 
and £o- Therefore, knowledge of B(z) in 1ow-k materials 
allows one to determine these two parameters. 

As predicted by Pippard and supported by the 
Bardccn-Cooper-Schricffcr (BCS) theory, in nonlocal su- 
perconductors B{z) is nonmonotonic with a sign reversal 
at a certain depth. This nonlocal electrodynamic effect 
is most pronounced in extreme type-I superconductors, 
such as In (k » 0.07); it can also occur in unconventional 
superconductors with nodes in the energy gap @ . 

The nonlocal effect has been directly confirmed by low- 
energy muon spin rotation (LE-^iSR) measurements on 
Pb, Ta and Nb 0,11], and soon after by measurements of 
polarized neutron reflectivity (PNR) in In Q. In Refs. @, 
8] a first attempt was made to infer £o from B{z) in Pb. 
However, due to an unresolved issue of systematic errors 
only statistical errors were estimated. One can resolved 
this issue combining LE-/^SR and PNR measurements 
since in the latter the systematic errors can be excluded. 

The knowledge of B(z) in nonlocal superconductors 
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also allows one to determine Z . Electrons near the Fermi 
surface are dressed by a cloud of virtual plionons, lead- 
ing to an enhancement of their effective mass and conse- 
quently to a reduction of the Fermi velocity v-p. Below the 
critical temperature T c the phonon mediated attraction 
of elections exceeding their screened Coulomb repulsion 
results in formation of the Cooper pairs. In BCS theory 
electron-electron coupling is weak: N(Q)V <C 1, where 
V is the pairing potential and -/V(O) is the single-spin 
electron density of states at the Fermi surface obtained 
from specific-heat measurements. Eliashberg's strong- 
coupling theory (SCT) is free from this limitation and 
agrees better with experiments [l|. 

In SCT the effective mass of electrons near the Fermi 
surface is m* = Zm and consequently the effective mass 
of the Cooper pairs is m* p ~ 2Zm. Correspondingly, Al 
and £o arc renormalizcd with respect to their values in 
the weak-coupling (wc) limit as 

A L = yfZ\t c , (1) 
io = GT/Z. (2) 

If Al(T — > 0) is measured, the factor Z can be inferred 
from Eq. ([1]). )^ c (0) can be obtained from the Faber- 
Pippard formula for A derived for the electrons not in- 
teracting with the lattice [l(|. This leads to 

where e is the electron charge. On the other hand, if 
£o is measured, Z can be calculated from Eq. (2) using 
the BCS definition ^ c = hv F /nA{0) = 0.18hv F /k B T c , 
where A(0) = A(T = 0) is the energy gap and h 
h x (27r) and k B are the Planck and Boltzmann constant, 
respectively. 

Unfortunately neither of these approaches is applicable 
for quantitative analysis since reliable calculation of v-p 
is not possible due to the complex topology of the Fermi 
surface in polyvalent metals PJJ. However, if both Al(0) 
and £o are known, one can eliminate Up, thus obtaining 

c 2 h 2 i 

12.57rT 2 e 2 7 ' A L (0) 2 £ 2 ' [ ' 

where it is taken into account that N(0) = 3'y/2ir 2 k B 
(7 is the electron heat capacity coefficient [H[). Eq. fl3J) 
is a general expression independent of the relationship 
between the current density and the vector potential and 
of the specific nature of the electron-electron pairing. 

In SCT Z = I + A m , where Am is a mass-enhancement 
parameter calculated in Ref. [l| using tunneling spec- 
troscopy experimental data. A m can also be obtained 
from McMillan's equation for T c [li^ . but this approach 
is less reliable [l|. Therefore, for nonlocal superconduc- 
tors Eq. (j4} furnishes a bridge between the B{z) and the 
tunneling data, hence providing an independent test for 
the validity of Al and £0 inferred from the B(z) data. On 



TABLE I: Parameters of the samples, d is the thickness, o is 
the oxide layer thickness (see Ref. ||), RRR is the residual 
resistivity ratio; I is the elastic mean free path. 



sample 


T C (K) 


size (cm) 


d (/jm) 


(nm) 


R.R.R 


I (fim) 
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FIG. 1: Time spectra of muons with energy of 22keV at an 
applied field of 47 Oe in the normal state (a) and in the Meiss- 
ner state (b). The curves are fits to the Gaussian (a) and the 
time-domain (b) models. The inserts present the field spectra 
P(B). In the insert in (b) bg is the background contribution, 
while the curve represents two Gaussian peaks to approximate 
the experimental spectrum. Ao is the maximum observable 
asymmetry of the muon decay. 



the other hand, the vast majority of superconductors are 
local. For those Eq. (|4]) allows one to determine £0 from 
the measured values of Al(0) and Z using, e.g., LE-/xSR 
and tunneling experiments, respectively. 

In this work B(z) is measured using the LE-/xSR [HI, 
HH and PNR [r| techniques. Combining these comple- 
mentary techniques enables an independent verification 
of the inferred values of £0 and Al (T) and an estimate of 
there total uncertainties. The samples were two indium 
films whose parameters are listed in table 1. The films 
were deposited on an oxidized silicon wafer (IN-1) and 
on a sapphire crystal (IN-2) held at room temperature 
via thermal evaporation of indium shots (99.9999% pu- 
rity) at a base pressure < 5- 10 -9 mbar. The IN-1 sample 
was used for PNR measurements and for LE-/iSR mea- 
surements down to 2.9 K. The IN-2 sample was used for 
LE-/iSR measurements at lower temperatures. 

The samples' DC magnetization exhibits a clear-cut 
first-order phase transition with deep supercooling. The 
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60 90 120 150 
depth (nm) 

FIG. 2: Reduced field B/Bo obtained from the LE-/iSR mea- 
surements using the Gaussian model. The depth is the aver- 
age stopping distance of the implanted muons calculated from 
the Monte Carlo TRIM.SP code. The solid (dashed) line is 
the field profile calculated from the nonlocal (local) theory at 
T = 2.92 K with £ =380 nm, A L (0)=30 nm and the dead layer 
20= 4 nm, inferred from the x 2 an d the time-domain analysis. 
The insert presents the same field profiles on a linear scale. 



measured critical field H C (T) and T e perfectly agree with 
values reported in literature [lj| [l7|. The surface of 
the films consists of nearly atomically flat terraces (root- 
mean-square roughness of the terrace areas is < 2 A) with 
a typical size around 5 /xm with voids in between. This 
size is much larger than £o and the total area of the voids 
does not exceed 3% of the sample surface. Therefore the 
terrace surface structure should not affect the electro- 
dynamic properties of the films. For both samples the 
elastic mean free path t >> £ j therefore the samples are 
type-I superconductors in the pure limit 0. 

The LE-/iSR experiments were performed at the 
/xE4 beamline of the Swiss Muon Source at the Paul 
Scherrer Institute fl8j . Typical spectra of muon polar- 
ization P are presented in Fig. 1. In the normal state B 
is uniform and the P{B) spectrum (insert in Fig. la) has 
Gaussian shape due to nuclear dipolc broadening. In the 
Meissner state B is not uniform and P(B) is not Gaus- 
sian. However, Gaussian spectra (curve in the insert of 
Fig. lb) provide a good first approximation. 

B(z) points obtained using a Gaussian model Q for 
the depolarization of the precessmg muons are presented 
in Fig. 2 along with the B(z) curve calculated from the 
local and nonlocal theories . The sample temperature 
was determined in situ based on the H C (T) phase di- 
agram obtained from magnetization measurements. As 
can be seen in Fig. 2, the "Gaussian" S-points exhibit 
a pronounced non-exponential depth dependence, consis- 
tent with the nonlocal effect. Discrepancies between the 
points and the "nonlocal" theoretical curve are mainly 
caused by incomplete adequacy of the Gaussian model. 
On the other hand, the qualitative consistency of the 
"Gaussian" points with the nonlocal theory justifies the 
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FIG. 3: x7 NDF 

versus £o and the corresponding polynomial 
fittings normalized to the minimum value for each temper- 
ature. The shaded area indicates the range of the best esti- 
mates for £o- The inserted table gives positions of the minima 
and their average with the standard deviation. 



application of a time-domain model fl9j , directly assum- 
ing the "nonlocal" shape of B{z) with Al as an adjustable 
parameter. 

In nonlocal superconductors the effective penetration 
depth A c ff ~ (A L £o)^ 0. Therefore either Al or £ has to 
be determined using an additional criterium. Since A c ft is 
mainly sensitive to Al and £o weakly depends on temper- 
ature, this parameter is £o- Its optimal value was found 
from x 2 analysis of the global fits (simultaneous fitting of 
all implantation energy data sets at each temperature). 
The fits were performed with different £o chosen around a 
theoretical value of 377 nm [H]. The graphs for x 2 / NDF 
(NDF stands for number of degrees of freedom) versus £o 
are presented in Fig. 3. Values of £o at the minima and 
the best estimate for £o arc given in the insert. The val- 
ues of Al obtained from the global fits with £o = 380 nm 
are presented in a summary graph in Fig. 5. 

The PNR measurements were performed on the 
D3 reflectomcter at the NRU reactor in Chalk River. 
The experimental data and simulations for the reflectiv- 
ity of neutrons polarized parallel (R + ) and antiparallcl 
(R~~) to the magnetic field and for the spin asymme- 
try (R + - R~)/(R + + R~) are presented in Fig. 4. In 
agreement with our LE-/xSR results and previous PNR 
results 0, the neutron spin asymmetry (Fig. 4, insert 
(a)) simulated with B(z) calculated from the nonlocal 
theory matches the experimental data significantly bet- 
ter than the simulation based on the London field profile. 
The "nonlocal" B(z) was calculated with £o = 380 nm 
obtained from the LE-/iSR data. The best match of the 
simulated spin asymmetry with the experimental data 
was achieved for A L (T = 0.3 K) = 28.0 ± 2.5 nm. 

Results for Al inferred from the LE-/zSR and PNR 
measurements are shown in Fig. 5. The values of Al(T) 
are consistent with each other and agree with the two- 
fluid formula 0. The best estimate of Al(0) is 30 ±2 nm. 

Having determined Al(0) and £o we calculated Z from 
Eq. Q . The results are summarized in the tabic inserted 
in Fig. 5, from which it can be concluded that the value of 
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FIG. 4: Neutron reflectivity in the Meissner state at T = 
0.3 K and Bo = 275 Oe. Insert (a): Measured spin asymmetry 
(points) and its simulation with B(z) obtained from the local 
(dashed curve) and nonlocal (solid curve) theories with Al = 
28 nm and £o = 380 nm. Insert (b): Measured (points) and 
simulated (curve) reflectivity in the normal state. 
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FIG. 5: Temperature dependence of the London penetra- 
tion depth in indium. The inserted table gives the values 
of Al(T — >• 0), £o and Z obtained in this work, and the values 
of Z calculated from the tunneling spectroscopy data [lj and 
from McMillan's equation for T c 



Z obtained in this work agrees with the value obtained 
from electron tunneling data. The latter implies that 
our results quantitatively confirm the nonlocal electrody- 
namic effect and confirm the validity of Eq. ((4]). 

In conclusion, high resolution measurements of the 
magnetic field profiles B(z,T) were performed, for the 
first time, in an extreme type-I superconductor (indium) 
combining LE-/1SR and PNR measurements. The B(z) 
profiles quantitatively confirm the Pippard/BCS nonlocal 
electrodynamic effect. The B{z) data were used to deter- 
mine the Pippard coherence length £n, the renormaliza- 
tion factor Z for the electron-phonon mass enhancement 
of the Cooper pairs, and the London penetration depth 
Al (T) . The experimentally verified equation ^ allows 
one to infer £o m any superconductor from the results 
for A L (0) and Z. 
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